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1 Introduction 

From the observational data (Type la supernovae [1], the cosmic microwave background 
(CMB) [2], galaxy distribution [3], and weak lensing [4]), we find that the ACDM model is 
the most generally accepted approach in cosmology. It is based on the Friedmann-Lemaitre 
solution of the Einstein equations with a cosmological constant A, which describes the uni- 
verse with a uniformly distributed dust. The observation also reveals the existence of large 
scale structures, i.e., the large voids with the scale of about lOOMpc as well as the clus- 
ters of galaxies. It is important to explain how such cosmological large scale structures are 
formed from the initially very small density fluctuations at the recombination epoch. The 
time evolution of fluctuations depends on the background universe, which is described by 
cosmological parameters such as the Hubble constant Hq, the density parameter r^o, and a 
cosmological constant A. Hence, the observations of the large scale structures as well as the 
CMB anisotropics give the constraints on the parameters. There are many works on the 
formation of large scale structures and their observational properties related to cosmolog- 
ical parameters. However, they are mostly concerned with over-densities, such as clusters 
of galaxies. Their dynamical characteristics act on their appearance in the redshift space. 
For example, because of the Doppler shift associated with the peculiar velocities of galaxies, 
the velocity dispersion within a cluster makes its shape elongated toward the observer, and 
hence it is named "Finger of God" [5] . The Kaiser effect is also caused by peculiar velocities 
and arises from coherent motions as the galaxies fall inwards towards the cluster center [6]. 
The Kaiser effect usually leads an apparent flattening of the structure, in accordance with 
the particular dynamics of the situation. 

On the other hand, the under-dense regions such as voids have not been intensively 
investigated as large scale structures as it has been done for over-dense regions. A few 
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models of voids have been investigated (see e.g. [7], and references therein) in order to 
interpret observations at high redshift. The hot and cold spots detected in the CMB may 
be associated with super-clusters and huge voids, which are identified in the Sloan Digital 
Sky Survey catalog [8]. Subsequently, the integrated Sachs- Wolfe effect [9] due to spherically 
symmetric super-structures suggested the presence of nonlinear voids or clusters with a radius 
100 ~ 200/i-iMpc [10]. 1 

Herein, we focus on the under-dense regions that undergo cosmological expansion with 
the aim to investigate the effect of the cosmological parameters on their dynamics. As a first 
step of our investigation for such under-dense regions, we discuss the case of a single isolated 
void embedded in the uniform and isotropic universe, which is the most simple example to 
understand. 

In §.2, we briefly introduce our model with the basic equations for a single void dynamics 
based on general relativity (GR). Our analysis on the expansion of the void is given in §.3, by 
which we confirm the previous result based on a covariant formulation of Newtonian dynamics 
[13]. In §.4, we analyze the shape of a void in the redshift space by which we show that such 
a distortion effect could be observable, similarly to "Fingers of God" or to the Kaiser effect. 
Conclusion follows in §.5. 

2 Modeling a spherical void in the expanding universe 

In this paper, we consider a single isolated void in GR to investigate its dynamics. We 
propose a model of a spherical void in the universe by using the thin-shell approximation 
based on the Israel's formalism [14]. A void in our model is described by an under-dense 
region filled by a uniform distribution of dust, which is embedded in the Friedmann-Lemai'tre 
(FL) universe with a positive cosmological constant A > 0. Here we will analyze only the 
case of an empty void. The extension to the case of a non-empty void is straightforward (see 
Appendix A). The equations of motion for a void were originally derived by Maeda and Sato 
[15], and later they were rewritten in a more convenient expression [16] that we will use here. 

2.1 Basic Equations 

We assume that a void is characterized by a spherical shell S of the radius R embedded in an 
expanding uniform background universe. It defines the boundary layer between two spatial 
regions filled by two different isotropic and uniform distributions of dust. For convenience in 
writing, we use the indices "— " and "+" for the variables related to the inner region and to 
the outer region, respectively. When no ambiguity is present, we will use the symbol "it" in 
order to write two equations as a single one. We assume that the corresponding space-times 
are described by a Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 

{sinxzt {k± = l) 
X± {k± = 0) (2.1) 
sinhxzt {k± = -1) 

^One interesting proposal for "dark energy" is an inhomogeneous universe model without any exotic ansatz 
such as introduction of dark energy or modification of gravity [If, f2]. In this model, we are living near the 
center of the under-dense region (a huge void). This model is fascinating if it turns out to be able to explain 
all observational constraints including the CMB anisotropy, although it is against the Copernican world view. 
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where k± are the signs of the scalar curvatures and a± the scale factors of the FLRW space- 
times. Both space-times satisfy the Friedmann equations 



-H± + ^ = ^^P±, with H± = — , and p± = pvac + P± (2.2) 
3 a± 

where the dotted variables denote their time derivatives, oc a^^ are the matter densi- 
ties, and pvac = A/(87rG) accounts for a cosmological constant (viewed as a contribution to 
gravity in term of energy density). These space-time regions are joined through S, the three- 
dimensional wor/d volume of S in the space-time {i.e., the hyper-surface that corresponds to 
the envelope of S evolving over time), which metric reads 

ds| = -dr^ + R^{T)dn^ (2.3) 

where r is the proper time of S. The continuity of the space-time requires the following 
conditions to hold 

fi = a+/+(xf ) =a-/-(x^''^) (2.4) 
dr^ = dtl - aliU) [dxf^y = dtl - al{t.) [dx^^^)" (2.5) 



These imply the following equations; 



^ = /± (X± j .± .± = a± — , / (x±) = yi - ^ (2-6) 

7+(/>+ + F+i?) = l-{f-v- + H^R) (2.7) 

7_ 1 , , 

' 7±= ; 2.8 

where v is the peculiar velocity^ of S and 7 the Lorentz factor. The jump condition on the 
extrinsic curvature of S provides us with the following basic equations; 

— — = --i+v+H+ H 2 (2-9) 

dt+ ^ ms 

= [47ri?VV"^]^ (2.10) 



^4j{f' + vHR)]^ = -^, (2.11) 

where ms = ms{t-\-) is the mass of the shell and the bracket [ ]^ stands for [^]^ = ^+ — 
(see [14-16] for details). 

The dynamics of the shell S is defined by Eqs. (2.6)-(2.10) with the constraint (2.11). 
In the case of an empty void (p^^ = 0), we do not have to solve Eqs. (2.7) and (2.8), which 
give the information for V- and t-. We use the constraint (2.11) to give an initial value of 
ms as well as to check the precision of the numerical integration schema. 



^It stands for the the radial component of the pecuhar velocity field, but we use such a shorter terminology 
thanks to its spherical symmetry. 
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2.2 Model parameters 

Since the outer region characterizes the background universe, the index "+"will be omitted, 
which enables us to retrieve the notations used in cosmology {e.g., the Hubble parameter of 
the universe H stands for the expansion rate of the outer region, and so on. 

By setting the present value of the scale factor = 1, the redshift of an event which has 
happened at time t and is observed today {t = to) is given hy z = l/a{t) — 1. The integration 
of the above differential equations requires initial values for the variables 

V, n- = —, S-^'-^-l, and A=^. (2.12) 

The first two (the radius of the shell R and its peculiar velocity v) describe the motion of 
S, the second two (the ratio of two expansion rates Ti' and the density contrast 6~) give 
the boundaries conditions between the inner and the outer regions, and the last two (the FL 
parameters and A) describe the background cosmological model. Herein, the initial values 
of the variables v and T-L~ are given at redshift z = Zi (i.e., at time t = ti). We choose them 
as follows : 

Zi = 100, Vi = 0, and 'Hi=l. (2.13) 

It turns out that the dependence of the expansion history of S on Vi and T-L^^ is negligible as 
long as these initial data correspond to a sufficiently early epoch^ (typically, Zi > 10). In order 
to compare the theoretical predictions with observations, the remaining model parameters 
are settled to the present values at redshift z = 0, i.e., Rq, 5q , Qq and Aq, by means of an 
iterative method^. 

In addition to the above parameters, we introduce the curvature parameter 

K = n + X-1, (2.14) 

which is used to classify the FL cosmological models. It stands for the dimensionless ex- 
pression of the scalar curvature of the comoving space, which is identified to an integration 
constant in the Newtonian dynamics [13]. 

Because the outer region of a void is described by the cosmological model, in order 
to analyze the effects of cosmological parameters on the proper dynamics of voids, we have 
first to study their contributions to the background cosmological dynamics. Although all 
cosmological parameters intervene to define the dynamics at large scale, it is known that 
the universe undergoes three phases; matter dominated(if 7^ 0), curvature dominated 
(if K 7^ 0), and vacuum-energy dominated (if A 7^ 0) expansion phases, depending on the 
relative magnitudes of cosmological parameters {e.g. see [17]). In the present investigation, 
we assume A > and hence the Hubble parameter H decreases with time and eventually 
toward its asymptotic value 

H^= lim H=J^. (2.15) 

a— >oo V o 

^Similarly as in the linear perturbation theory, once the "decaying mode" becomes negligibly small, the 
eventual dynamics does not depend on initial conditions. 

''Namely, the initial values Qi and Xi are adjusted to find the present observed values (e.g., = 0.3 and 
Ao = 0.7) by solving the Friedmann equation for the outer region. We also solve the Friedmann equation for 
the inner region, which gives the relation between the variables of the inner region and those of the outer 
region by assuming the value 5(7 . Then, we solve the evolution equation for R with time. The initial radius 
Ri is also adjusted by means of an iterative method to obtain the present size of a void Rq. 
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In that case, there are two distinct behaviors which are characterized by the sign of k ( see 
Eq. (2.14)): Either H decreases monotonically toward Hoo (if k < 0) or it reaches first a 
minimum -ffmin < -f^oo and it eventuaUy grows toward Hoc from downward (if k > ) . This 
minimum value is found to be 



Hmin — Hoc\ 1 — . ° 9 ^ (2.16) 

at a = l^lo'^o ^) where kq is the present value of the dimensionless curvature constant. Near 
this minimum, the universe experiences a loitering period. 

3 Dynamics of voids 

With initial conditions given in Eq. (2.13), we suspect that, a spherical void tends to grow 
because the under-dense inner region expands faster than the outer region.^ To describe such 
a process in the Newtonian dynamics, instead of the peculiar velocity v, we have used the 
dimensionless quantity y = R/HR, which enables us to disentangle/define a "proper motion" 
of the shell S from the background expansion velocity. Here in order to confirm the result in 
the Newtonian dynamics, we discuss the same variable y in GR, which reads 

_ dR/dt f'v 

y=^ = '^HR- ^'-'^ 

The variable (y — 1) works as the correction factor of the peculiar velocity v of the shell S 
to the Hubble expansion velocity in the outer region. 

The proper motion of a void can be found from Eq. (3.1) in the comoving space, which 
frame is adapted to the investigation of large scale structures since one obtains a genuine 
interpretation of their dynamics^. With this schema in mind, we investigate the dynamics of 
an empty spherical region in the universe. As for the non-empty case, we shall discuss it in 
Appendix A. 

We analyze the behavior of a spherical empty void embedded in the background universe 
characterized by fio = 0.3 and Aq G {0,0.7,1.4}, which represents distinct cosmological 
models. For each one, the evolution of y with time has been calculated for two spherical 
voids. Their radii at initial time, which we set Zi = 100, are chosen to be 10% and to 100% 
of theinitial horizon size; i.e., Ri = wH^^ with w G {0.1, 1}. The results are described in 
Fig. 1, where the correction factor y is given in terms of the redshift (or the scale factor a). 
We investigate the effect of A on the void dynamics by comparing the Aq = 0.7 case (which 
will be used as a standard) with the other ones. Our analysis, which refiects the standpoint 
of the observer since it is based on the measurable quantities, is as follows : 

1. The voids always grow faster than Hubble expansion (i.e., y > 1) along their evolution. 
The peculiar velocity starts with a huge burst and decreases asymptotically toward 
Hubble velocity; the present value is higher than the Hubble flow by about 10%. 



^In the Newtonian dynamics of an empty spherical region, it is understood by the fact that each point of 
S is radiaUy accelerated by the attractive force of its closest denser outer region, because of the symmetry of 
its shape. 

^Indeed, when the co-existence of the CMB blackbody spectrum with the recession of galaxies is para- 
doxical with respect to Newtonian mechanics, these observations are reconciled within a description of a 
thermodynamical equilibrium in GR [18]. The set of points in the space-time from where the CMB can be 
observed with a black body spectrum at a given temperature To defines the comoving space and provides us 
with a universal reference frame as long as FL model is in agreement with observations. 
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2. The smaller the radius gets, the higher the peculiar velocity is. Namely, the void with 
the largest initial size Ri = H~ grows slightly slower, which provides us with the 
minimal peculiar velocity. 

3. While the dependence on Aq is not significant in the beginning, its effect eventually 
appears in the evolution, and the correction factor y increases with Aq. With Aq = 0.7, 
the peculiar velocity reaches the value higher than the Hubble expansion at z ~ 10 by 
about 20%. 

4. The redshift z* when y reaches its maximum r/max = y{z*) decreases with Aq. It 
corresponds to z* ~ 40 for Aq = and to z* ~ 1.7 for Aq = 1.4. For the intermediate 
value Aq = 0.7, for which the corresponding curve shows a plateau, A gives the maximal 
contribution to v at ~ 1.7, which corresponds to 30% of that of the matter density 
(when it is compared to the case of Aq = 0). 

5. Let us pay attention to the presence of a bump on the curve, which becomes visible 
for Aq > 0.7 {i.e. when > 0). It is caused by the fact that the universe experiences 
a loitering period of cosmological expansion, while the void then continues its own 
expansion. 

6. The present values yo related to = 0.3 but to different values for Aq G {0, 0.7, 1.4} 
are very close to each other, while the variation of y with time (or with z) depends 
undeniably on Aq. In other words, the A effect, which accounts for a deviation between 
these curves, is not substantial nowadays. 

7. It is remarkable that the present GR approach for the sub-horizon void confirms the 
previous result found in the Newtonian dynamics [13]. Even for a relativistic spherical 
void with a horizon size radius Ri = , the relativistic effect turns out to be weak. 

To disentangle a A effect on the dynamics of a void is quite a tricky task since all the 
model parameters intervene in its dynamics. It is a matter of fact that the voids expansion is 
magnified with Aq for a fixed value of Q,q. At first glance, such a feature could be understood 
as a gravitational local effect similar to the one described by the Schwarzschild solution of 
the Einstein equations with a cosmological constant A. Actually, this interpretation is an 
artifact due to the comparison of the A effect between different cosmological models defined 
by the curvature parameter kq. The A effect intervenes on the expansion of the outer region 
at large scale and not as being part of the own void dynamics. Such a property can be 
understood from a dimensional analysis of the Einstein equations, which shows that the 
presence of A in the gravitational field is felt only at great (cosmological) distances [19]. As 
a result, the excess of the peculiar velocity v of an empty void with respect to the Hubble 
expansion velocity depends mainly on fio [16], which is related to the density of the outer 
region. Such a behavior is similar in the case of an under-dense region but weighted by the 
density contrast 5q . With respect to the observer's standpoint, the dependence on Aq, which 
intervenes on the Hubble expansion, remains extremely weak at high redshift z > 20 and low 
redshift z ~ 0, but it appears prominent at redshift z ~ 1.7. At this redshift, the peculiar 
velocity v for Aq = 0.7 (i.e., kq = 0) is magnified by up to 50% compared with the case 
of Aq = {i.e., kq = —0.7) (see Figl). Such features characterize Aq through the behavior 
of the cosmological expansion H. If A > then H decreases with time and will reach its 
asymptotic value Hqo- There may be a loitering period during which a spherical void can 
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Figure 1. The correction factor to the Hubble flow — For two voids with initial radii Ri = Q.1H~^ 
(solid line) and (dashed line), the correction factor y is depicted in terms of the redshift of 

the background universe z by assuming three cosmological models defined by fio — 0.3 and Aq G 
{0,0.7,1.4}. 



grow much faster compared to the cosmological expansion. It clearly appears as a bump in 
the curve for Aq = 1.4 in Figl. Such a feature could be used to evaluate the value of Aoby 
a statistical analysis of cosmological parameters from observational data [13]. Here we will 
show one simple example, by which the peculiar velocity of a void shell creates an observable 
deformation of a void shape in the redshift space. 

4 Redshift of a void 

To analyze the large scale structures in the universe, one makes use of the redshift- space. It is 
identified to the direct product of the observer celestial sphere and the redshift of sources z. 
However, the images of structures are deformed because the redshift, although it is used as 
distance indicator, does not provides us with an exact value of the distance. As a result, the 
image of a spherical void located at cosmic distance does not appear spherical in the redshift 
space'^. The shape of a void reflects its dynamical characteristics. Namely, the redshift of 
photons emitted from galaxies at the edge of the void (the shell S) takes into account its own 
expansion as well as the evolutionary history of the universe.^ 

Hereafter, we restrict our analysis to a spatially flat spaces of the FLRW world model, 
which is characterized by a single parameter fig = 1 ~ -^O) in addition to the Hubble constant 

'^The redshift space is often used by observers because it is one of the simplest ways to construct the large 
scale structure at cosmological distance scale from observational data. 

®A sphere which expands in the comoving space with the cosmological expansion (but not a static one) 
provides us with a deformed image in the redshift space because of the evolutional effect of the universe. As 
we shall see later, the peculiar velocity of a void shell enhances this tendency and its deformed shape includes 
the information about cosmological parameters. 
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Hq. In Sec. 4.1, the shape of a void in the redshift-space is modeled. We focus to an empty 
void {p^^^ =0), in which the deformation reaches its maximal magnitude (see Appendix B). 
A qualitative analysis is performed by numerical calculations in Sec. 4.3. 

4.1 Redshift by proper motion of a shell 

The void appears as the intersection of the observer's past light cone with the world volume 
S of the void shell S. We shall visualize such an image of the intersection by means of 
the redshift of the source on S and its angular position on the celestial sphere. With the 
Euclidian spatial coordinates, the FLRW metric reads 

ds^ = -dt^ + a^{t){dx^ + dy'^ + dz'^) . (4.1) 

Prom the symmetry, without loss of generality, it is enough to analyze the light rays propa- 

y 




Figure 2. The light rays from a void shell to the observer. — The observer O is located at the 
origin (0,0), the x-axis stands for the line of sight toward the void center (xv,0). Its boundary layer 
(the shell S) is represented by the circle of comoving radius r = r{t). The light ray emitted from 
a point {xs,ys) located on the backside of S with respect to the observer at comoving distance /, 
propagates inside the void before crossing S* at a point x , and then reaches the observer O; ip and 
a are the viewing angles of the light source with respect to the x-axis from the observer O and from 
the center of S, respectively. The light ray from the foreside point {x'^, y'^) propagates directly toward 
the observer O. 



gating on a two-dimensional plane, which we choose the xy-plane (see Fig. 2). The observer 
O is located at the origin of the coordinates system and the center of the void shell lies at 
X = xy on the x-axis. The comoving radius of the void is given by r(t) = fix^^\t))i the 
peculiar expansion velocity is u = adr/dt, and the Lorentz factor is 7 = l/\/l — v"^. A point 
XgS on S is described by its spatial (comoving) coordinates as Xs = {xs,ys,0), where 

Xs = Xy + r cos ip, ys = rsinip, (4.2) 

and the 4-velocity of S at Xg reads 

7, — cos if, — sin 93, . (4.3) 

a a / 
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The 4-momentum of a light ray reads 



^ / cos a sma \ , , 

F = E 1, , ,0 . 4.4 

\ a ay 

where E is the photon energy and a is the angle between the light ray and the x-axis (see 
Fig. 2). The 4-velocity of the observer reads Uo = (1, 0, 0, 0). The redshift of a photon 
emitted from Xg and received by the observer is then given by 

l + zs = = 1^7 [1 + cos{^ - «)] (4.5) 

As we shall see later, the geometrical effects inside the void are negligible. Hence, instead of 
solving the geodesic equation to calculate Eg/ Eg, we use the (homogeneous) approximation 



Eg _ Eg 
Eo Eo 



FLRW '^s 



l + z, (4.6) 



where z stands for the (cosmological) redshift when the photon is emitted. This is free from 
the Doppler effect due to the shell expansion. Eq. (4.5) transforms 

l + Zg = {l + z)-i[l + vcos{ip-a)] (4.7) 

The comoving coordinates of the emission event read 

dz 

Xg = l{z)cosa, T/s = /(z) sina, l{z) = / — (4-8) 



aoH{z) 

where l{z) is the comoving distance from the observer to Xg. Hence, for a given a, from 
Eqs. (4.2) and (4.8), we find the redshift z when the photon is emitted and two intersecting 
points (two values of (p that correspond to points located in the foreside and in the backside 
of S with respect to the observer). Since v is derived from the dynamics equations, we obtain 
Zg in terms of a, which relation provides us with the shape of a void in the redshift space. 
Its numerical procedure is summarized as follows: 

(i) Solve our basic equation for the void dynamics to find r = r(z) and v = v(z). 

(ii) For a given angles a, find the intersecting points {f) of the light path with the shell 
and the time z when the photon is emitted. 

(iii) Calculate the redshift Zg of the void shell for each angle a. 

(iv) Draw the shape, i.e. {zx,Zy) = {zg{a) cosa, Zs{a) sin a). 

The order of magnitude of the contribution to redshift Zg from the peculiar velocity v 
is estimated by the related Doppler shift as^ 

^•^slooppier ~ cos((^ — a) 0.1ffi?cos((/9 — a), for r2 = 0.3. (4-9) 

^It must be emphasized that the interpretation of the redshift z due to the Friedmann expansion as the 
Doppler shift is quite tricky without a model which enables us to disentangle the Doppler effect by the motion 
of an object from the gravitational redshift. 
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4.2 Other effects on the redshift of a void 

In the previous subsection, we have evaluated the redshift of a void just due to the Doppler 
effect by the proper motion of the sheh. However, there are two other effects which contribute 
to the redshift. Those are the following two effects caused by density inhomogeneities. 



• The Sachs- Wolfe effect. — The Sachs- Wolfe effect vanishes [20] because there is no 
potential difference between the shell S and the background universe under the thin- 
shell approximation. The non- vanishing redshift is found for the integrated Sachs- Wolfe 
(ISW) effect [21], which is evaluated by 

Az.lisw ~ {HRf (4.10) 

• A gravitational lens effect. — For the light ray from the backside of a void, the light 
is bended at some point on the foreside shell (the mark 'x" in Fig. 2). The scattering 
angle at the shell (see Appendix B) is approximately given by 

Aa^0.25HRsm{(p-a), for JIq = 0.3 (4.11) 

Note that Aa = for a photon which goes through the center ((/? = a = 0) of the void 
or for one which comes from the edge ((/? — a = 7r/2). Hence the ratio of the maximum 
radius to the minimum one is not modified, although this light bending may slightly 
change the shape of a void in the redshift space. 

In the present thin shell void model, the only relevant quantity is related to the change 
of the optical path length inside the void caused by a deflection, which is ~ i?(Aa)^. 
Hence, the redshift due to the gravitational lens effect reads 

R(Aa)^ r 

Azs\gl^— ^-0.06 — {HRfsm{(p-a), for fio = 0.3 (4.12) 

As a result, the contribution to redshift of these inhomogeneity effects turned out to be 
negligible compared to that of the Doppler effect given in Eq. (4.9) for voids with sub- horizon 
sizes. 

4.3 Numerical results 

For the analysis of a spherical void in the redshift space, we compare the image of its shell S 
with the one of a comoving sphere of constant comoving radius r, which we use as a reference 
shape. In the following figures, the image of S is represented by the intersection of S in the 
redshift space with the light cone of the observer. The event in the redshift space is described 
by the coordinates (z cos a, z sin a) in the two-dimensional xy plane. As we shall see from 
our numerical results, the image of S (red dots) looks roughly like an ellipse whose center is 
located at the redshift zy. Its major axis (of length 2Az||) coincides with the line of sight of 
the observer, and the minor axis (of length 2A2:_|_) is free from the elongation. As a result, 
the shape of a void in the redshift space is elongated in the direction of the line of sight. 
A parameter to quantify such an effect is the deformation ratio^^ 

n = . 4.13) 



^°It is the opposite deformation of Kaiser effect, since tlie peculiar velocities are oriented outward. 

^^In the present case, it is given by the ratio of the semi-major axis to the semi-minor axis of the approximate 
ellipse since A^n > Azx. It is more appropriated than the eccentricity for measuring the elongation of the 
image because it shows whether S is expanding (JL > 1) or contracting (7?. < 1). 
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The shape of the reference void is also an ellipse (black dots), which is due to the evolutionary 
effect of the background universe. It shares the minor axis with the image of S because there 
is no redshit due to the Doppler effect in that direction, and shows a smaller deformation 
ratio (i.e., it is less elongated). The comparison between these images provides us with an 
estimate on the magnitude of the peculiar velocity of the void shell (relatively to the Hubble 
flow). 
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Figure 3. Empty spherical voids in the redshift space. — The images (red dots) of their boundary 
layers (i.e., the void shell S) are depicted for the case of = 0.3 and Aq = 0.7. The present values 
of their radii are (a) Rq = O.IHq^ and (b) Rq = O.OIHq^. The observed radii at zy = 1 are : 
(a) Ry = 0.0437iJo"^ and (b) Ry = 0.00437ifo-^ The images of the standard static spheres (small 
black dots) are also displayed. 



In addition to the above main results, we analyze the void shapes with respect to the 
following observational viewpoints : 

(1) Dependence of properties of voids 

• Size — For the case of = 0.3, we have analyzed the images of voids located at = 1 
with two different present radii Ro G {OAH^^, O.OIHq^},^^ which are depicted in Fig. 3. 
The image of the void with the present size of Rq = O.IHq^ shows a small fore-back 
asymmetry (see Fig. 3 (a)). Such a feature can be interpreted as an evolutionary effect 
due to the photons emitted at different times from points of the expanding shell S. It 
does not substantially affect the shape. It is not visible on the image of the void with 
a smaller radius Rq ~ O.OIHq^, for which it is expected to be weaker because the time 
difference is small (see Fig. 3 (b)). Note that these shapes are similar while the radii of 
corresponding voids differ by one order of magnitude. Therefore, the deformation ratio 
TZ turns out to be (almost) independent of the size of a void, for which the elongation 

The observed void size is 

Rvizv) = aor{zv)^Rox'^^ <Ro, (4.14) 

r(0) 

because the photons were emitted at z = zy (the past). 
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in the redshift space is proportional to Az<j ~ ~ O.IHR. Such an interesting feature 
should be taken into account in the observation of voids. 



Distance — For f^o = 0.3, the images of voids with the resent radius Rq = O.IHq ^ but 
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Figure 4. The images of empty spherical voids with the radius Rq = O.IHq in the redshift space 
located at two different redshifts (a) zy = 0.5 and (b) zy — 2. The observed radii of voids are (a) 
Ry = 0.0616Hq^ and (b) Ry = 0.0269i/(^^, respectively. The small black dots represent the images 
of the standard static spheres 



located at two different redshift zy G {0.5,2} are depicted in Fig. 4. We find that the 
deformation ratio TZ increases with redshift. Such an effect is clearly seen in Fig. 5 (a), 
in which four images of voids with Rq = O.IHq^ at redshift zy G {0.5,1,1.5,2} are 
displayed. The reason why this elongation of the shape with the redshift (distance) zy 
increases is twofold, i.e., in the redshift space, Az± gets smaller because it corresponds 
to the observed size Ry{zy), while A^n increases with the peculiar velocity of the void 
shell (see footnote 11). 

(2) Dependence of the background cosmological models 

• Matter density — Four images of voids with Rq = O.IHq^ located at redshift zy G 
{0.5, 1, 1.5, 2} are displayed in Fig. 5 for two different values of Qq G {0.3, 1}. 

By comparing these images in Fig. 5 (a,b) , it is clear that the deformation ratio TZ 
increases with i}Q for all redshift. Let us remind that the background matter den- 
sity, which characterizes the void as well as the outer space where it is embedded, is 
represented through the cosmological parameter Qq. 

• Deformation ratio — The property of the deformation ratio TZ of being independent on 
the void size motivates us to analyze this parameter with more care. In Fig. 6(a), we 
confirm that Az± is a monotonically decreasing function of zy. On the other hand, Az|| 
increases for = 1 while it first decreases toward a minimum and eventually increases 
in the case of fio = 0.3. The deformation ratio TZ is an monotonically increasing function 
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Figure 5. The images in the redshift space of an empty spherical void with present radius Rq 
O.IHq^ at redshift zy e {0.5, 1, 1.5, 2} for (a) = 0.3 and (h)no = 1. 



of zv because and hence the peculiar velocity of the void shell {v ~ 0.2HR), increases 
with zy. This tendency is enhanced as l^o increases (see Fig. 6 (b)). Since we focus on 
a flat universe model (k = 0), any information for Q provides us with an estimate of a 
cosmological constant. 

5 Conclusion and Remarks 

With the aim to understand the dynamics of voids in the universe and the effect of a cos- 
mological constant A on it, we analyze a relativistic model of a spherical under-dense region 
in a FL universe. The inner and the outer regions are described by space-times with uni- 
formly distributed dust. We study the effects of cosmological parameters on the dynamics 
of voids, focusing on the case of an empty void. We have calculated the peculiar velocity 
of a void shell, with which the shell moves in the comoving space because the lower density 
region expands faster than the outside Hubble flow. We show that the dynamics depends 
mainly on the density parameter Qq. Our analysis confirms the previous result obtained in 



- 13 - 



0.-15 

0.4 
0.35 

0.3 
0J5 

02 
0.15 

0.1 
0j05 





1 1 1 1 


^,1 








^11. 


- 

















i7 — rt 11 




—9- m--—^ 


















- + -+j 


-+ - ; 




1 1 1 1 




"■'"■•="-wrr , 



0.5 



1.5 



(a) Az_L and Azu 




(b) The deformation ratio TZ 



Figure 6. (a)The semi-major axis Az|| and tlie semi-minor axis Azj^ in terms of the distance zy. 
We find that Az± always decreases as zy increases because the observed void size gets small. On 
the other hand, Az|| increases as zy increases due to the Doppler effect except for nearby voids for 
f^o = 0.3. (b) The deformation ratio TZ of the shape of a void in the red shift space, which is defined 
by 7?. = Az||/A2:_L, where Az|| and Az± are the semi-major axis and semi-minor axis radii of the 
void shape in the redshift space, respectively. The ratio TZ increases as the distance zy increases. We 
also find the ratio TZ highly depends on the background cosmological model, which may give us some 
information on the cosmological constant. 



the Newtonian dynamics, namely, the relativistic effects are rather weak. The shape of voids 
in the redshift space has also been investigated in the case of a spatially flat cosmological 
model. It is characterized mainly by the void kinematics, i.e., the Doppler shift due to the 
proper motion of a void shell. The other effects on the shape, such as the ISW effect and the 
gravitational lens effect, turn out to be rather small for a void with a sub-horizon size. The 
void in the redshift space appears as an ellipse shape elongated in the direction of the line 
of sight (opposite to a deformation by the Kaiser effect), although a small modification from 
the ellipse shape by an evolutionary effect may become perceptible with increasing the size. 

Since the void dynamics depends on a background cosmological model, such observa- 
tional data could be used to improve statistics of cosmological parameters, especially that of 
a cosmological constant by means of the structures at redshift z ~ 2. It is, however, clear 
that our model of a single spherical void is too simple and it must be improved in order 
to account for shape irregularities due to interactions of galaxies on the edges of voids. In 
order to perform such improvements, one may need numerical approach. Since the dynamics 
of such a void-network is much more complex than that of a single isolated void, A^-body 
simulation would be required. In the conventional A^-body simulation, we usually assume 
that the space-time is described by a uniform and isotropic universe. The expansion velocity 
of the background universe is uniform and it depends on the mean density of matter fluid but 
does not on the local density. Such an approach is justified if the density perturbations are 
small. However, if the perturbations are non-linear, i.e., if the density contrast is very large 
just as in the late stage of structure formation, we may have to worry about the position 



- 14 - 



dependent background expansion just as our present void model. In fact, according to our 
model, which describes the dynamics of a spherical under-dense region in a FL universe, such 
an under-dense region expands faster than the background Hubble flow, i.e., the space-time 
region with under-dense matter distribution expands faster than the over-dense region. This 
fact will change the dynamics of structures just as our simple case. Hence, with the aim to 
understand the mechanism of the large scale structure formation, we have to improve the 
Newtonian A^-body simulation. For this purpose, the Newton- Friedmann model [13], which 
is based on the covariant formulation of Euler-Poisson equations [22], appears to be more 
adapted for being implemented in the numerical codes, although we may need further study 
in its practical formulation. 

It is interesting to note that such an investigation can be extended to test other gravity 
theories [23] or approaches for dark energy. It is widely discussed as an alternative interpreta- 
tion of the observed acceleration of the cosmological expansion [1-4], although a cosmological 
constant remains the most appropriated explanation according to Occam's razor. 
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A Spherical under-dense void 

In this appendix, we extend our analysis of the void dynamics to the case of non-empty 
void. We assume that the inner region of a void is filled by a uniform distribution of dust 
with a density p^™\< p*-™^). The density contrast is characterized by 6~ = p^^Vp*'™'' ~ 1 
(— 1 < 5~ < 0). We restrict our analysis to the case of a flat background universe (k = 0). 
With appropriate initial values for the model parameters, we have integrated the evolution 
equations for an under-dense void, which variables are (y,Q,6~), until the redshift zq. If 
zo = 0, the calculation gives the present values. It turns out, however, that the relation 
between our variables {y,i},d~) does not depend on zq as well as the initial time Zi as long 
as the "decaying mode" becomes negligibly small during the evolution. 

The relation found by the numerical calculation is shown by circles in Fig. 7. We can 
find a fitting formula for the relation as 

y-1 = — = — +0.1|<5-|2 + 0.07|<5-|=^), for n + X = l (A.l) 
HK 

which is shown by the solid curved in Fig. 7. It approximates our numerical results very well. 
The difference between the numerical results and the formula (A.l) is less than one percent. 

The expansion of an under-dense region differs from that of an empty void {6~ = — 1), 
although the dynamics is similar to that of an empty void. As we see from Fig. 7, for a 
given 0, the higher the inner density is (i.e., 5~ — )• 0), the smaller the peculiar velocity of 
the void shell (i.e., (y — 1)) is. This is because the space with under-dense matter expands 
faster than that with over-dense matter. 

As we have seen in Fig. 1, the peculiar velocity increases monotonically as f2 increases. 
The similar feature is seen for non-empty void (see Fig. 7). In fact the curves in Fig. 7 are 
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Velocity y versus and 5 . 

Figure 7. The velocity y of a non-empty void. — We assume + A = 1 and 6^ E {01, —0.5, —0.1}. 
The circles and the continuous lines correspond respectively to our numerical results and to the fitting 
formula given in Eq. (3.2). 



found by extrapolating the curve in Fig. 1 with (reasonable) continuous transformations up 
to the constant function y = 1 for 6~ = 0. As the result, we find that the correction factor 
y depends strongly on 6~ as well as on il.. 

B Deflection of the backside light ray at the boundary shefl 

In this appendix, we derive the deflection angle for a light ray that is emitted from a point 
located on the backside of a void with respect to the observer. As described in Fig. 2, it 
propagates inside the void, crosses the boundary shell at a point denoted x, and then 
reaches the observer O. For convenience, we perform the coordinate transformations of the 
4-momentum in the vicinity of the shell four times as follows : 

1. We translate Eq. (4.4) from Cartesian coordinates (4.1) to spherical ones (2.1), i.e., 

E E 

k^ = E, k^ = cos((f-a), /c^ = — sin(v9-a), /c^ = , (B.l) 

a af 

2. We adopt a Gaussian normal coordinate system (r, n, 9, ip) which covers both the outer 
and the inner parts of S. n is the spatial coordinate normal to the hypersurface of the 
shell S, which corresponds to n = 0, and the 3-metric of S is given by Eq. (2.3). The 
coordinates transformation of A;'^ at n = from the outer FLRW frame to the Gaussian 
normal frame gives 

E 

k'^ = E^[\+v cos{ip — a)], k^ = —E^[v+cos{ip—a)], /c'^ = — — sin((^ — a), /c^ = , 

R, 

(B.2) 
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where we have used the relation (see [23]), 

dt dr dt dr 'J 

OT OT a on on a 

3. The transformation of from the Gaussian normal frame to the inner FLRW frame 
gives 

/c*" = £'7+7_[l + (u+ — V-) cos((/? — a) — f+W-], 

= [— cos((y9 — a) — V-^- + + v^v- cos((/7 — a)l, 

a_ 

A;'^ = — ^sin((/5-a), A;'' = . (B.4) 

This gives the boundary conditions for null geodesic equations in the inner region. 

In general, a local scattering angle is not well defined when the spatial curvatures of the 
outer and the inner parts are different. However, in the case of an empty void (p^^ = 0), 
the inner region is described by a de Sitter space-time and the spatial curvature is negligible. 
Hence, for this case, we can introduce another Cartesian coordinates system 

ds2 ^ _^^2 _^ a'i{t-) {dxi + dy'i + dz^_) , (B.5) 

and perform a coordinates transformation of from the spherical coordinates to the Carte- 
sian coordinates as 

E 

k^~ = — [— cos a + {v- — f +) cos (/?] -|- 0(vj.) 
E 

F- = —[-sma + {v. - v+) sin v?] 0(u±) . (B.6) 
Finally we obtain the direction angle and the scattering angle as follows 



^ / \ sm(99 — a) ^, 2\ 

tana_ = = tana -|- (f+ — v^) ^ h U[v<) 

k^- cos^ a 

Aa = a_ — a = (f+ — v^) sm{Lp — a) -|- Oivj^) . (B.7) 
According to our numerical analysis, we find as follows: 

• If 17o = 1, then v+ ^ 0.2HR and ^ -0.3HR . 

• If 17o = 0.3, then v+ ^ O.IHR and V- ^ -0.15HR . 

Hence, we obtain the upper bound of the scattering angle, which is given in Eq. (4.11). 
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